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We report first principles calculations, within density functional theory, of copper intercalated 
titanium diselenides, Cua;TiSe2, for values of x ranging from to 0.11. The effect of intercalation on 
the energy bands and density of states of the host material is studied in order to better understand 
the cause of the superconductivity that was recently observed in these structures. We find that 
charge transfer from the copper atoms to the metal dichalcogenide host layers causes a gradual 
reduction in the number of holes in the otherwise semi-metallic pristine TiSe2, thus supressing the 
charge density wave transition at low temperatures, and a corresponding increase in the density of 
states at the Fermi level. These effects are probably what drive the superconducting transition in 
the intercalated systems. 



I. INTRODUCTION 

Transition metal dichalcogenides (TMDCs) are quasi- 
two-dimensional, highly anisotropic layered compounds 
that are of great interest in basic research as well as in 
applications in such areas as lubrication, catalysis, elec- 
trochemical photocells, and novel battery systems^. Re- 
cently, it became possible to synthe size other forms of 
TMDCs such as balls and nanotubes^lil. TMDCs gen- 
erally have the formula MX2, where M is a transition 
metal, such as vanadium, titanium, tantalum, molyb- 
denum, or others, and X is a chalcogen atom, such as 
sulfur, selenium, or tellurium. Each layer consists of a 
metal sheet sandwitched between two chalcogen sheets. 
Whereas within each layer the bonds connecting a metal 
atom in one sheet to the chalcogen atoms in surrounding 
sheets are strong covalent bonds, adjacent layers are only 
weakly coupled through the van der Waal's interaction. 
Two types of MX2 sandwitches occur depending on the 
coordination of the transition metal atom by the chalco- 
gens. In IT-MX2 the coordination is octahedral whereas 
in 2H-MX2 it is trigonal prismatic. 

The weak coupling between adjacent layers in TMDCs 
makes it possible to introduce foreign species into the 
region between these layers, a region called the van der 
Waal's gap. This process, known as intercalation, pro- 
vides the ability to tune the electronic properties of the 
system by controlling the type and concentration of the 
intercalant. This is a general property of layered ma- 
terials, and just like graphite, which has been interca- 
lated with many kinds of atoms and molecules, 5 TMDCs 
have been intercalated with alkali atoms, transition metal 
atoms, and organic and inorganic molecules^. 

Among the most important properties of TMDCs is 
that some of them undergo a phase transition to a charge 
density wave (CDW) state as the temperature is low- 
erecP. lT-TiSe 2 is one of the first compounds wher e the 
CDW state was observed at a temperature of 200K 6 7 8 . 
The state is associated with the formation of a commen- 
surate (2x2x2) superlattice. It is not completely clear 
what mechanism drives this phase transition. Possible 



explanations are offered in terms of an indirect Jahn- 
TellereffectpESl or by an exciton condensation mecha- 
nism™ 

Recent angle resolved photocmission spectroscopy 
(ARPES) studies have not resolved whether lT-TiSe 2 
is a semimetal with low electron and hole concentraion 
or a narrow-ga p semic onductor in either the normal or 
the CDW stateP^SEIl Optical measurements, on the 
other hand, reveal that TiSe 2 is a semimetal with low car- 
rier density at temperatures above as well as below the 
CDW transition^. This is consistent with band struc- 
ture calculations 1 which reveal the presence of an elec- 
tron pocket, around the L point in the Brillouin zone, 
derived from Ti 3d states, and a hole pocket around the 
r point derived from Se 4p states. 

Recently it was reported^ that upon controlled inter- 
calation of TiSc 2 with Cu to yield Cu^TiSea the CDW 
transition is continually suppressed, making way for a 
new superconducting state that emerges near x = 0.04, 
with a maximum transition temperature T c of 4.15K 
ocurring at x = 0.08; further increase in x results in 
a drop of T c . It should be noted that while superconduc- 
tivity was previously observed in 2H structures, such as 
2H-NbSe 2 , 2H-TaSe 2 , 2H-NbS 2 , and 2H-TaS 2 ™the cop- 
per intercalated TiSe 2 is the first superconducting 1T- 
structured TMDC. As Cu atoms are intercalated, they 
move into the van der Waal's gap separating adjacent 
TiSe 2 layers. With a single electron in the outer 4s shell, 
electronic charge transfer proceeds from Cu to the TiSe 2 
layers, causing an increase in the number of electrons in 
the Ti 3d bands and a reduction in the number of holes in 
the Se 4p bands. As a result of this, the effective screen- 
ing of the electron-hole Coulumb interaction is enhanced, 
causing a drop in the CDW transition temperature as the 
intercalant concentration increases, and at x — 0.06 the 
transition is completelly supressed. The superconduct- 
ing transition, on the other hand, emerges at x = 0.04. 
The CuzTiSe 2 phase diagram resembles the correspond- 
ing diagram in cuprate and heavy fermion superconduc- 
tors, except that there it is the antiferromagnetic spin 
order that competes with superconductivity, whereas in 
Cu x TiSe 2 it is the charge order^Sl. 
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Recently, calculations of density of states (DOS) and 
the electron phonon coupling, using linear muffin-tin or- 
bital (LMTO) method within the local density approx- 
imation, were reported for the case x — 0.08 by Jeong 
and Jarlborgpl The authors consider a 2x2x3 supercell, 
Cu(TiSe2)i2, and conclude that charge transfer from cop- 
per raises the electronic density of states at the Fermi 
energy. They also find a strong dependence on pressure 
of the DOS, and consequently the superconducting tran- 
sition temperature. 

In this paper we report band structure calculations 
for Cua;TiSe2 for values of x ranging from 0.0 to 0.11 
in an attempt to study the effect of copper intercalation 
on the electronic structure of TiSe 2 . In particular, we 
calculate the energy bands of TiSe2 and the compounds 
Cu 1 / 16 TiSc 2 , Cu 1 / 13 TiSe 2 , and Cu 1 / 9 TiSe 2 . In Section 
we discuss the computational method and some relevan 
details, and in Section III we discuss the results and con- 
clusions. 



II. METHOD 




FIG. 1: (color online)The unit cell of TiSe2 is shown in (a), 
while the Brillouin zone, and the k-path along which the en- 
ergy bands are calculated, are shown in (b). In (c), the unit 
cell for Cu(TiSe2)g is shown. 



The first-principles calculations presented in this work 
were performed using the all-electron full potential lin- 
ear augmented plane wave plus local orbitals (FP- 
LAPW+lo) method as implemented in WEN2K code^l. 
In this method the core states are treated in a fully rela- 
tivistic way but the valence states are treated at a scalar 
relativistic level. The exchange-correlation potential was 
calculated using the generalized gradient approximation 
(GGA) as proposed by Pedrew, Burke, and Ernzerhof 23 . 

In the calculations on lT-TiSe 2 the experimentally 
measured lattice constants (a = b = 3.534A, c = 6.008A) 
are employed^. The distance betwen the Ti plane and 
either of the adjacent Se planes is d — zc, where the 
value of z is obtained by minimizing the forces on the 
atoms to below 2.0 mRy/Bohr. In calculations on the 
intercalated compounds Cuj;TiSe 2 the lattice constants 
are taken to coincide with the experimental values^, but 
the internal structure of the unit cell is allowed to relax, 
subject only to the constraint that the space group re- 
mains unchanged; again the force on each atom is reduced 
to below 2.0 mRy/Bohr. For the compound Cui/i 6 TiSc 2 
the structure is taken to be that of TiSe 2 4x4 superlat- 
tice in the a-b plane, with one Cu atom placed midway 
between one Ti atom in one layer and the nearest Ti 
atom in the adjacent layer; the unit cell is then hexago- 
nal with a' = b 1 = 4a. Accordingly, the unit cell has 1 
Cu atom, 16 Ti atoms, and 32 Se atoms, and if one Ti 
atom is placed at the origin (0,0,0), then the Cu atom 
is placed at (0,0,c'/2), where c' is the lattice parameter, 
in the direction normal to the a-b plane, in the interca- 
lated systems. A similar situation holds for Cu!/9TiSe 2 
and Cu 1 / 13 TiSe 2 where the corresponding superlattices 
in the a-b plane are 3x3 and V4~3xV'i3, respectively. The 
space group is P-3ml(#164), the same as in TiSe 2 , in 
all these structures, except in Cu 1 / 13 TiSe 2 , where it is 



P-3(#147). It should be noted here that the superlattice 
structure employed in these calculations is not necessar- 
ily the experimental lattice structure. For example, for 
the case x — 0.0625, we assume that the system is pe- 
riodic with lattice constants in the a and b directions 
that are 4 times bigger than the corresponding values 
for x — 0, whereas the experimental crystal structure 
is not reported and there may well be some disorder in 
the system such that the Cu atoms, residing between ad- 
jacent TMDC layers, do not form a perfectly periodic 
two-dimensional lattice. 

For calculations in this work, the radii of the muffin- 
tin spheres for Cu, Ti, and Se atoms were taken to be 
2.5ao, 2.5ao, and 2.25ao respectively, where ao is the Bohr 
radius. The number of k-points in the irreducible part 
of the Brillouin zone (BZ) were 38, 41, 47, and 360 for 
Cu 1 / 16 TiSc 2 , Cu 1 / 13 TiSe 2 , Cu 1 / 9 TiSe2, and TiSe 2 , re- 
spectively; these numbers were chosen so that the density 
of k-points in the whole BZ is roughly the same in the 
four cases. For all structures considered in this work we 
set the parameter RmtKmax=7, where Rmt is the small- 
est muffin-tin radius and Kmax is a cutoff wave-vector. 
The valence wavefunctions inside the muffin-tin spheres 
are expanded in terms of spherical harmonics up to l max 
= 10 while in the interstitial region they are expanded 
in plane waves with a wave- vector cutoff Kmax , and the 
charge density is Fourier expanded up to Gmax— 14a^ 1 . 
Convergence of the self consistent field calculations is at- 
tained with a total energy convergence tolerance of 0.1 
mRy. 



FIG. 2: (color online)Band character plot for TiSe2 along high symmetry directions in the Brillouin zone. In (a), the contribution 
from Se 4p states is shown. The width of the band lines is proportional to the contribution of the Se 4p states to these bands. 
In (b), the contribution of the Ti 3d states to the energy bands is emphasized. The Fermi level is at zero energy. 
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TABLE I: Comparison of the crystal structure of Cu(TiSe2)g 
before optimization, as derived from the stucture of TiSe2, 
with the relaxed structure obtained by minimizing the forces 
on the atoms. The space group of the crystal is P-3ml. 



III. RESULTS AND DISCUSSION 

It is well known that the calculation method followed in 
this work (GGA, PBE) overestimates the unit cell volume 
and lattice constants and it gives a very large c/a value in 
layered compounds 25 . Indeed, upon optimization, Titov 



et alP^ obtained a value for c/a in TiSe2 which is about 
5% larger than the experimental value. Since the lattice 
constants are known for the systems we study here, we 
perform the calculations using these experimental values, 
but we allow the internal coordinates to relax. 

The unit cell of TiSe2, shown in Fig.[lR, consists of one 
Ti atom and two Se atoms. The crystal coordinates of 
these atoms are given by Ti(0, 0, 0), Sel(l/3, 2/3, z), and 
Se2(2/3, 1/3, —z). Upon minimizing the forces on each 
atom to below 2 mRy/Bohr, it is found that z=0.2573; 
the distance between the Ti plane and an adjacent Se 
plane is thus 1.546A. The Brillouin zone (BZ) corre- 
sponding to the TiSe2 structure is shown in Fig. [TJd, and 
the calculated energy bands along some high symmetry 
directions in the BZ are shown in Fig. [2j In Fig. [2]i the 
contribution to the bands from the Se 4p states is shown, 
while in Fig. [2|i the contribution from the Ti 3d states is 
displayed. The bands are in broad agreement with previ- 
ous calculations^. They are also in good agreement with 
recent ARPES study^ 3 which shows that lT-TiSe 2 has a 
scmimetallic band structure with a negative energy gap 
Eq = -70±15 meV, the Se 4p band having a maximum 
at 30±10 meV and the Ti 3d band having a minimum at 
-40±5 meV with respect to the Fermi energy. By exam- 
ining these plots, it is observed that the electron pocket 
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FIG. 3: (color online)Energy bands for the 4x4x1 TiSe2 supperlattice are shown in (a), while those for Cu(TiSe2)i6 are shown 
in (b). The contribution of Ti 3d states is emphasized in both plots, the width of the band lines being proportional to this 
contribution. The Fermi level is at zero energy. 



at L is derived mainly from Ti 3d states, while the hole 
pockets at T and A are derived from Se 4p states with 
some significant hybridization with the Ti 3d states. 

Upon intercalation by copper, the c-axis is increased 
slightly, and some slight rearrangement of the positions 
of the Se atoms takes place. The new positions are found 
by relaxing the structure so as to minimize the forces act- 
ing on the atoms. The lattice constants are held fixed at 
the experimental values, while the internal coordinates 
are allowed to change. We find that the incorporation of 
copper atoms between TiSe2 layers leads to only a slight 
modification in the atomic positions. To illustrate this 
point, we present in Table [T] the coordinates of the in- 
equivalent atoms for the case of Cu(TiSe2)g in the ideal 
superlattice structure (before optimization) and the cor- 
responding coordinates after optimization. In the three 
compounds, Cu x TiSe2, for x = 1/16, 1/13, and 1/9, the 
Se atoms are displaced in the z-direction by Sz where 
-0.006A < Sz < 0.009A, whereas they are displaced in 
the a-b plane by less than 0.037A. Because of the slight 
displacement in the z-dircction of the Se atoms, they no 
longer all lie in a single plane as is the case in undoped 
TiSe 2 . Nevertheless, our results indicate that the dis- 
turbance caused by intercalating copper into the van der 



Waal's gaps of TiSe2 is small as far as its effect on the 
atomic positions is concerned. This is not unusual, for 
a similar situation holds in graphite intercalation com- 
pounds where, besides an increase in the c-lattice con- 
stant, the effect on the structure of the carbon planes is 
minimal. This is a consequence of the fact that, both in 
TiSe 2 and graphite, the bonding within each layer ( in 
TiSe 2 a layer consists of three atomic sheets) is of the 
strong covalent type. 

The intercalation of copper into the van der Waal's 
gap in TiSe2 is stabilized by charge transfer from cop- 
per to the TMDC layers. Here we calculate the charge 
on each atom using Bader's quantum theory of "atoms 
in molecules" (AIM)PS1. In this theory space is divided 
among the atoms so that each atom occupies a definite 
volume bounded by a surface at which the flux in the 
gradient of the electron density p(r) vanishes, 



h(r) ■ Vjo(r) = 



where n(r) is a unit vector normal to this bounding sur- 
face at r. The number of electrons in a given atom, oc- 




FIG. 4: (color online) Energy bands for the \/T3x\/T3xl TiSe2 superlattice obtained by zone folding of the TiSe2 bands, and the 
energy bands for Cu(TiSe2)i3, are shown in (a) and (b), respectively. The Fermi level is at zero energy. The Ti 3d contribution 
to each band is proportional to the width of the band line. Note that while the Fermi level rises by only 0.06 eV, as a result of 
Cu intercalation, the hole band, with a maximum at 0.17 eV at T in (a), shifts down in energy to just below Fermi level. 



cupying a volume fi, is then given by 

N e = [ p(r)d 3 r. 

For undoped TiSe2 we find that Ti is positively charged 
with a charge of 1.575e while each Se has a nega- 
tive charge of -0.787e. In the intercalated compounds 
Cu(TiSe2) x we find that the copper atom acquires a pos- 
itive charge of 0.445e, 0.444e, and 0.450e, for x = 9, 13, 
and 16, respectively. To illustrate how the charge is dis- 
tributed we present in Table |TT] the charge on each atom 
in the case of Cu(TiSe2)g, as calculated by Bader's AIM 
theory. 

For the copper intercalated systems, the energy bands 
appear to be far more complicated than in undoped TiSe2 
because the unit cell has many more atoms. A direct 
comparison between the undoped TiSe2 energy bands 
and those of the doped compounds is not very helpful 
in clearly elucidating the effect of copper intercalation. 
Most of the bands in the intercalated compounds result 
from the zone folding, into the smaller BZ, of the bands 
in the larger BZ of the pristine compound. Therefore, in 
order to gauge accurately the effect of intercalation with 



copper on the energy bands, it is more useful to compare 
the bands in the Cu intercalated TiSe2 with those that 
result from zone folding the pristine TiSe2 bands into the 
smaller BZ characteristic of the intercalated compounds. 
For example, the compound Cui/igTiSe2 has the chemi- 
cal formula Cu(TiSe 2 )i6, and except for a slight variation 
in the c lattice parameter and an extra Cu atom, its unit 



Atom Site Q/e AQ/e 

Ti(l) la 1.551 -0.024 

Ti(2) 6g 1.563 -0.012 

Ti(3) 2d 1.577 +0.002 

Se(l) 6i -0.820 -0.033 

Se(2) 6i -0.801 -0.014 

Se(3) 6i -0.796 -0.009 

Cu lb 0.445 +0.445 



TABLE II: The charge on each inequivalent site in 
Cu(TiSe2)g, as indicated by AIM theory, and the correspond- 
ing change, resulting from copper intercalation, in the charge 
on each site. Note that the charge of an electron is -e, so that 
e is positive. 
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FIG. 5: (color online)Energy bands for the 3x3x1 TiSe2 superlattice, obtained by zone folding of the TiSe2 bands, are shown in 
(a), while the bands of the copper intercalated compound is shown in (b). The Fermi level is at zero energy. The contribution 
of the Ti 3d states is emphasized. Note the rise, by 0.08 eV, of the Fermi level in (b), as compared to (a). The hole band, with 
a maximum at 0.17 eV at V in (a), shifts down in energy to -0.08 eV, as a result of copper intercalation. A similar fate meets 
the hole band at A. 



cell would be the same as that of a 4x4x1 superlattice 
of pristine TiSe 2 , which also consists of 16 TiSe 2 units. 
To understand the effect of the addition of copper in this 
case, we compare the energy bands of Cu(TiSe2)i6 with 
those of (TiSe2)i6, the 4x4x1 undoped TiSe 2 superlat- 
tice. Similar arguments apply in the case of intercalation 
compounds for other values of x, the copper concentra- 
tion. 

The energy bands of the systems studied in this work 
are shown in Figs. [3] [4] and [5j In Fig. [3] we compare 
the energy bands of the 4x4x1 superlattice with those of 
Cu(TiSe2)i6- It should be noted that for the 4x4x1 super- 
lattice, the BZ is 16 times smaller than that for pristine 
TiSe 2 . The M and L points in the pristine TiSe 2 BZ are 
zone- folded into the T and A points, respectively, in the 
superlattice BZ; this explains why features in the band 
structure at the M and L points in Fig. [2] now appear 
at points T and A, respectively, in Fig. [3] For example, 
in Fig. [3^, the band along TA, which crosses the Fermi 
level, and which is derived from Ti 3d states, corresponds 
to the band along ML in the pristine TiSe 2 BZ; the M 
point is folded to T and the L point is folded to A. The 
effect of intercalation with Cu is seen to be a rise in the 



Fermi level due to charge transfer from Cu. The bands 
derived from Cu states are more than 1 eV below the 
Fermi level; thus the main influence of Cu intercalation 
on the electronic properties is a charge transfer to the 
host layers, since in the immediate vicinity of the Fermi 
surface, all the bands are derived from Ti and Se states. 
For the case of the 3x3x1 and vT3x%/l3xl super lattices, 
the M and L points in the bigger BZ, corresponding to 
the case when the unit cell contains only one TiSe 2 unit, 
are folded, respectively, into the M an L points in the 
smaller BZ of the superlattices. 

By closely examinning the energy bands in Figs. [3]-|5] 
we conclude that the charge transfer from Cu to the host 
layers leads to a slight rise in the Fermi level, but the host 
bands derived from Ti 3d states are largely unaffected; 
the rigid band model is a very good approxiamtion in 
describing the effect of Cu intercalation in TiSe 2 as far 
as the bands derived from Ti 3d states are concerned. 
For the Se 4p derived hole pockets at T and A points, 
however, the effect is not described by the rigid band 
model. The slight displacement of the Se atoms as a 
result of intercalation, and the charge transfer from Cu, 
combine to cause a suppression of the hole band. For 
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FIG. 6: Electronic density of states for TiSe2 and two interca- 
lation compounds. Note the increase in the value of the DOS 
at the Fermi level, set at the zero of energy, as a result of Cu 
intercalation. 



example, for copper concentration x—1/9, we see from 
the energy bands in Fig. [5] that while the Fermi level 
is raised by ~0.08eV due to Cu intercalation, the hole 
pockets at T and A are lowered by ~0.25eV; indeed, these 
hole pockets are now below the Fermi level, and thus they 
are completely filled. Within the exciton condensation 
mechanis m 1 11 ! 12 ! for CDW transition, such a reduction in 
the number of holes leads to a suppression of the CDW 
state. 

The emergence of superconductivity, as the concentra- 
tion of Cu increases, is probably related to the suppres- 
sion of the CDW transition and the concomitant increase 
in the electronic density of states in the vicinity of the 
Fermi level. We calculated the density of states for the 
various systems considered and in Fig. [6] we plot the den- 
sity of state for pristine TiSe2 and two of the three inter- 
calation compounds considered in this work. According 
to the BCS theory of superconductivity,^ the transition 
temperature T c increases with an increase in the density 
of states in the immediate vicinity of the Fermi level, the 
other important factor affecting T c being the strenght of 
the electron-phonon coupling. In Table |III[ we present 
the calculated density of state at the Fermi energy, in 
states per eV per one TiSe2 units, for pristine TiSe2 as 
well as the other various intercalation compounds. Note 
the increase in the density of states resulting form inter- 



calation by Cu; the calculated result that this density is 
largest for the Cui/i3TiSe2 case is consistent with the ex- 
perimental observation that the superconducting transi- 
tion temperature T c is largest for Cu concentration given 
by x w 0.08. However, as pointed out earlier, the density 
of states is not the only factor affecting the value of the 
transition temperature. 

It should be remarked that the density of states at 
the Fermi energy, N(E^), is not necessarily an increas- 
ing function of x, the copper concentration, as Table III 
clearly indicates. The value of N(E^) is intimately con- 
nected to the band structure. As x increases, the Fermi 



Compount 



N(Ef 



TiSe 2 

Cu 1/16 TiSe 2 
Cui /13 TiSe 2 
Cu 1/9 TiSe 2 



0.86 
1.06 
1.22 
1.01 



TABLE III: The calculated electronic density of state at the 
Fermi level N(Ef), in the states per eV per one TiSe 2 unit, 
for TiSe 2 and its copper intercalation compounds. 



level rises; consequently, the number of bands crossing 
the Fermi level may change, or that a band will cross the 
Fermi level at a new K-point in the BZ where DOS may 
have a larger or a smaller value than the old K-point. 



IV. CONCLUSION 

In conclusion, we have presented calculations, within 
the density functional theory, of the electronic energy 
bands and the density of state of copper intercalated 
TiSe 2 . We have shown that a rigid band model describes 
reasonably well the effect of the intercalation on Ti 3d de- 
rived bands, but not as well for Se 4p derived bands. As a 
consequence of charge transfer from the Cu atoms to the 
host layers, the Fermi level rises, the Se 4p hole pocket 
gets filled; this is the probable cause for the suppression 
of the charge density wave transition. This, coupled with 
an increase of the electronic density of states at the Fermi 
level, is perhaps the cause of the emerging superconduc- 
tivity in these novel intercalation compounds. 
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